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Abstract 

We propose a stochastic process for stock movements that, with 
just one source of Brownian noise, has an instantaneous volatihty that 
rises from a type of statistical feedback across many time scales. This 
results in a stationary non-Gaussian process which captures many 
features observed in time series of real stock returns. These include 
volatility clustering, a kurtosis which decreases slowly over time to- 
gether with a close to log-normal distribution of instantaneous volatil- 
ity. We calculate the rate of decay of volatility-volatility correlations, 
which depends on the strength of the memory in the system and fits 
well to empirical observations. 

1 Introduction 

In the past years there has been a body of work investigating the statistical 
properties of financial data [1, 2], in particular of time series of stock returns. 
Such data is available on the smallest of time scales, namely tick by tick, 
ranging to years and years. Certain properties appear to be extremely stable 
across different stocks and different time horizons. For example, histograms 
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of stock returns exhibit fat tails such that the cumulative distribution decays 
as a power law with exponent —3. This behaviour is extremely stable, decay- 
ing slowly towards a Gaussian distribution as the time scale of the returns 
is increased. Time series of returns exhibit bursts of correlated volatility, 
in the sense that there will be regions of higher or lower volatility, together 
with a persistence in the sense that the volatility exhibits a rather significant 
correlation over time. In addition, the distribution of empirical volatilities 
can be fit very well by a close-to log-normal distribution. 

The standard log-normal model of stock returns, inspired by earlier work 
due to Bachclicr [3] and used for example by Black, Scholcs and Merton [4, 5] 
for the purpose of option pricing, has had great success as the basic building 
block for modeling stock returns. However, none of the interesting features 
of financial time series which we mentioned above are captured by such a 
model. Instead, a variety of other models have been proposed in order to 
more realistically describe financial data [6, 7, 8, 9, 10, 11, 12]. For example, 
stochastic volatility models [6, 13, 14, 15] capture well the volatility cluster- 
ing and fat tails, except that these tails decay way too quickly towards a 
Gaussian distribution. Multifractal models have also enjoyed great success 
in describing many of the stylized facts [7, 8, 9], although they contain an 
additional source of noise, in contrast to the model we present in this pa- 
per. It is thus still an exciting open question, to probe the dynamics giving 
rise to price formation both on a fundamental level, namely by studying the 
dynamics of the order books [16, 17, 18] as well as on a more phenomenologi- 
cal level, namely by understanding the possible types of stochastic dynamics 
that could underlie price formation. The current paper will address this issue 
along the latter lines. 

We recently proposed modeling the driving noise of stocks by a statistical 
feedback process of the type presented in [19]. Based on that model, we were 
able to quite well capture many features of stock option markets, lending 
weight to the model [21, 20, 22]. However, as pointed out in those papers, 
the model is not an entirely realistic one. The main reason is that there is 
one single characteristic time in that model, and in particular the effective 
volatility at each time is related to the conditional probability of observing an 
outcome of the process at time t given what was observed at time t — 0. This 
is a shortcoming of that model for one of real stock returns; in real markets, 
traders drive the price of the stock based on their own trading horizon. But 
there are traders who react to each tick the stock makes, ranging to those 
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reacting to what they beheve is relevant on the horizon of a year or more, 
and of course, there is the entire spectrum in-between. Therefore, an optimal 
model of real price movements should attempt to capture this existence of 
multi-time scales. 



2 The Model of Stock Returns 
2.1 The statistical feedback model 

First we summarize the original model, in a slightly different notation. If S 
is the stock price, then y = InS" is the log-stock price. We proposed ([20]) 
that 

dy = iidt + adfl (1) 

with 

dn = p{n\noy-^duj (2) 

Here, /i is the instantaneous log return, cr is the volatility parameter, and 
a; is a standard Brownian noise such that {du}{t)du}{t')) = 6{t — t')dt. In 
all that follows we set = for simphcity. P is the probability density 
of fl conditioned on Qq (which can be chosen to be 0). Following the lines 
of standard stochastic calculus, P evolves according to a nonlinear Fokker- 
Planck equation of order 2 — q, namely 

BP 1 8^ 

dt 2dQ? ^ ' 

Exact time-dependent solutions of this equation can be found [23] and are of 
the form of Tsallis distributions of entropic index q (equivalent to Student 
distributions) 

P(Qi I Qo) = {a,V + ^(Q* - Qo)')^ (4) 



with 



t 



{{2-q)i3-q)c,r (5) 
q-1 



(2-q){3-q) 



(6) 
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and 7 = 1^^. The g-dependent constant is given by Cq — 



TT '■9-1 2' 

1-^ r2(-i^) 



From equation Eq (1) with // = it is evident that 

^t^- (7) 
a 

so in the following we refer to y as our variable of interest. 

The index q corresponds to the entropic Tsallis index within the general 
thermo statistics framework [24] . It is also related to the degrees of freedom 
in a Student distribution. If g = 1, we recover a standard log-normal process 
for stock returns. If g > 1, the effective volatility of the process is controlled 
by P raised to a negative power. Therefore, fluctuations are enhanced if 
extreme values of y are realized, and they are more moderate if y is less 
extreme relative to the information available at time 0. This gives rise to 
fat tails and bursts in the effective volatility of the process. Indeed, for 
q > 1 P oi Eq(4) exhibits heavy tails. In particular, g = 1.5 well-models 
the distribution of empirically observed log stock returns, reproducing the 
observed cubic decay of the cumulative distribution of returns over short 
time scales. 

The model here implies that the ensemble distribution of yt — yo is of the 
Tsallis form for all times. It is harder to say something about the distribution 
of increments over time yt+i — yt- However, if we plot out histograms of 
these increments we see that they have statistical properties consistent with 
empirical observations of log stock returns. But because the process explicitly 
depends on the initial value yQ at time 0, and therefore also explicitly on the 
time t, it is clearly not an optimal model of real returns, although it captures 
several realistic features. 

Written in a slightly different notation so that the time ti is denoted by 
i and the time at a discrete increment At later is denoted by ti + At — ti+i, 



Vi+i ^yi + (TiAui 



(8) 




(9) 
(10) 



where we have replaced yi = Qi/a and define 

aifl = aqa'^i"' 



(11) 
(12) 
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We denote 

Pgii, 0) = (a,,o + k^oiy, - yo)')^ (13) 

But for a multiplicative pre-factor, this is the conditional probability density 
of observing yi at time i given yo at time 0. In the notation above it is 
clear that the process with volatility a is simply a time- and state-dependent 
standard stochastic equation. It just so happens that the solution yields a 
probability distribution for the variable y that is of the Student or Tsallis 
form. 



2.2 The multi-time scale model 

As a more realistic approach, we propose the following model, a simple gen- 
eralization of the one we previously proposed. We write the volatility as 

j=i-N 

where N is the size of some past moving window which contributes to the 
effective volatility of the price, and 

= p.ihjy-' (15) 



\2 



= [a. 



+ bij{yi-yjr) (16) 



This leads to 



yi+i — yi 



+ ^( E -LO^A... (17) 



N j=i-N 



Here, Pq{i,j) is of the form Eq (13), with the argument i, replaced by i,j 
and where 

ai,j = a^a'^ii-jy (18) 
hi,j = \{t~3)-^ (19) 

Since market participants are reacting to information on many different 
time horizons, the effective volatility of the process for y is modeled here as 
the sum of contributions from all the different time horizons. In other words, 
from each time scale i — j, there is a collective feedback into the system 
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based on how extreme the movement of yi is perceived on the relevant time 
scale. If each of the feedback processes on the individual time scales i — j 
were independent, in the sense that only the dynamics across that time scale 
would feed back into the system, then Pg would be the probability density 
of the variable yi conditioned on an initial value yj at time j. (Note that 
one possible more general formulation of this model would be to replace 
Pq{i,j) with the true probability density P{yi \ yj) of the model. But it is 
non-trivial to solve for this quantity based on the non-Markovian situation 
at hand, and so we choose instead to study the current simplified scenario). 
Thus, across each time horizon, the contribution to the volatility follows the 
same dynamics as an independent statistical feedback process on that time 
scale. Our initial model of [21, 20] is recovered if only one time horizon is 
assumed relevant. 

In the current form, the memory in the system is described by an equal 
weighted moving average that does not decay. We can easily modify this, by 
introducing other types of averaging kernels Wij resulting in 

dyi = — ^ , .^^ cL}i (20) 

One possible choice of weights which we study in this paper is exponential 
decay Wij = exp(— A(i — j))- If we let ^ oo this model represents pure 
exponential decay. For finite N it describes a model with a cut-off exponential 
decay. 

Note that in this process there is only one driving noise acting at time i, 
namely the Gaussian random variable cjj. Via the statistical feedback pro- 
cess on different time-horizons, the effective volatility at time i is composed 
as the sum of a scries of apparently random variables each the result of a 
statistical feedback process on the corresponding time scale. Since we know 
the dynamics of each of these variables, the hope is that we can perform 
useful calculations related to the statistics of y. 
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3 Theoretical Calculations 



3.1 The volatility 

We shall see that it is straightforward to calculate the second moment, or 
volatihty, of this process. We insert a as in Eq(15) into Eq (17) and take the 
continuous time limit where 

i = ti-^t (21) 

j = tj^S 

N = tN^T 

to get 

1 /■* 1 
dyt = j^_^^^t,s + Ksiyt - ysf)ds)^duji (22) 

Integrating over all past times s which contribute to the process yields an 
effective volatility. We see that it must contain some sort of memory or 
volatility correlation. Explicitly, for the short time volatility correlation, 

{dytdyt^) = h{f {t-syaqa^ + hq{t-s)-\yt-ysf)dsY^ 

1 Jt-T 

( r ((^' - ^'sV^i'^'' + M-t^' - ^T\yt' - ys'f)dt'f^dujtdujt>) 

Jt'-T 

(23) 

where we have used the property that {dutdut') = 6{t — t')dt. 

We evaluate this integral self-consistently by replacing (y^ — y^)^ with its 
expectation. We obtain 

which does not depend on absolute time. The fact that squared increments 
are constants proportional to dt implies a normally diffusive stationary 
process. (Observe also that if one so wishes, can always be chosen as 
^22^-7 order to cancel out the T-dependent part of {dyl).) Our result is 
exact, and valid as long as 6, < 1. 
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Very similar results are obtained if we have an exponential decay rather 
than a sharp cut-off, as in Eq (40). In the continuous time limit keeping the 
notation of Eq(35) one obtains 

2\ _ (//_Texp(-A(t - g)(a^agP + V^Hl/t - ysf)ds)dt 
Jt-T exp(-A(t - s))ds 

Replacing {yt — inside the integral by its expectation, we solve self- 
consistently to obtain 

= ^y(7 + ^) dt = A,dt (26) 
where F is the standard Gamma function and T ^ oo. 



3.2 Volatility- volatility correlations 

We can also look at the correlation of squared price differences of this process, 
which in discrete notation can be written as 

{^yj/^yli) - (Ay5(Ay^+,) = ?,^t\{ajal,) - {aj){al,)) (27) 

If / = the two are perfectly correlated. If / = the two are perfectly 
uncorrelated (if we restrict our analysis to the accuracy of the leading order 
in g — 1, which means that we ignore the fact that i/i+iv itself depends on past 
values of y in proportion to a — 1 factor) . This is assumed to be the case 
in the discussion that follows (see also [25]). In that setting, the correlation 
decays from the one extreme to the other as the lag / increases from to N . 
Correlation is due to the fact that there is a common set of histories which 
both ^yf and ^yfj^i depend on. The quantity we are interested in is the 
autocorrelation of volatihties, namely Cj^i+j = {'^t^'i+i) ~ {^t) {^f-^i) ^ ^ seen 
here: 

i-l i-N+l i-1 i+l 

j=i+l—N j=i-N j=i+l-N j=i 

= {iC + D)iC + D))-{a^){al,) (28) 

where the notation should be obvious. D and D correspond to the contribu- 
tions of the volatility- volatility correlation which are uncorrelated with each 
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other and with both C terms. The C terms on the other hand are correlated 
with each other. 

In Eq(28), the terms a approach constants as — > 1, so the dependency 
on the random variables Ui — yj disappears and the problem becomes trivial. 
For g > 1, as / increases the common set of past values of yj which the two 
volatilities depend on decreases. This leads to a decrease in the volatility- 
volatility correlation at the rate U{1). In our case, when / = there is 
full correlation and this quantity is at a maximum. When I > N, there is 
complete de-correlation and this quantity should be zero. 

From Eq(28) follows that the uncorrelated contribution to {(yffyf^i) can 
be written as {Daf + DC). This is the amount that gets subtracted away 
from the ((Xj) {(Ji+i) term of the Cj_j_|_; equation. Therefore, the rate at which 
the volatility autocorrelation decreases can be expressed as 

U{l)^{ai){al^)-{ba^ + DC) (29) 

This expression can be evaluated explicitly both for constant and exponen- 
tially decaying weights, as presented in the Appendix. In Figure 4 and the 
following paragraphs, the results of this calculation, namely Eq (35) and Eq 
(40) together with Eq (36) Eq (39) are discussed in more detail. 

4 Simulations and Numerical Results 

We created numerical simulations of the process (always with cr = 1 for 
simplicity), an example of which is shown in Figure 1. From this plot, non- 
Gaussian features such as bursts of volatility are apparent. Not only do we 
see this volatility clustering by eye but also in the plot of Figure 2 where the 

variogram of volatility is depicted, defined as f = ((^ dy^ — dy^^^ ). For a 
random walk with no memory in the volatility, one would obtain a straight 
line. On the other hand, there is significant correlation in the volatility for 
q — l.h and an exponential decay. Qualitatively similar results are obtained 
for other g > 1, and other cut-off lengths, with constant or exponentially 
decaying weights. 

Another feature of stock returns is that the kurtosis approaches zero 
slowly as the time lag of the returns increase. Over short time-scales, re- 
turns have a high kurtosis (are very non-Gaussian with fat tails) but as the 
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Figure 1: A typical path of the multi timc-scalc non-Gaussian model with q ~ 1.5, equal 
weighted window with cut-off at N = 2000 (which corresponds to T = 140 days using 
dt = 1./14 as the simulation step) is shown. Very similar results are obtained for a wide 
range of cut-off lengths, and choice of exponential weight. 

time scale increases, the kurtosis tends to vanish and returns appear more 
Gaussian. The rate of this decay of kurtosis can roughly be approximated 
by U{1), the rate of decay of Ci^i+i [1]. In Figure 3 we plot the empirical 
kurtosis calculated from simulations of the process using an equal weighted 
window, together with a plot of U of Eq(35). Also shown is the hue depicting 
which we use as a proxy for empirical decay of the kurtosis of real stock 
returns [1]. It is clear that our model indeed has a very slowly decaying 
kurtosis. In fact, with an equal weighted window it decays even slower than 
what is observed for real stock returns (for lags sufficiently smaller than the 
cut-off time). However, if we include an exponential weight rather than a 
constant one, then the decay of the kurtosis is more realistic. 

In fact, it is easy to see that it is the rate of decay of the memory across 
time-scales which controls the rate of decay of the kurtosis. This is explored 
in Figure 4. The parameter which we vary is the half-time of the exponential 
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Figure 2: A volatility variogram shows the correlation of volatility across time scales. 
Here for q — 1.5 and an exponentially decaying weight with half-time at / = bOOdt « 36 
days and T — bOOOdt w 1 year. The solid line corresponds to a theoretical curve assuming 
the variogram goes as Eq(30) with g = l^'^-^'^ which is a proxy for the behaviour of real 
stock returns. The inset depicts the q ^ 1 case (standard Black-Scholes log- normal model), 
that shows no correlation or memory. 



weight defined as Zq.s = log{2)/\. For numerical reasons, we kept large yet 
finite in all simulations, at a value of = 5000, corresponding to T ^ 1 year 
with our choice of dt. Consequently, it is convenient to rescale I by T, and 
study the behaviour of the process as h = I0.5/T varies from 00 (a constant 
equal weighted scenario) to (no memory at all). The rate of decay of the 
kurtosis varies from that of Eq(35) which is close to constant for times less 
than the cut-off time T, to roughly which is the result we expect for 
perfectly uncorrelated Gaussian random variables. From the plot in Figure 
4b it seems that choosing a value of h in the range 0.02 = 0.1 reproduces a 
decay of kurtosis reasonable when compared to that of real stock returns, if 
we look at the region / -C T. These values correspond to half-times of 7 days 
and 36 days, respectively, with the choice of simulation parameters used. 
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Figure 3: The average kurtosis calculated from (top curve) a simulated time series of 
the process with constant equal weights, together with the U{1) Eq (35). The kurtosis 
decays slower that that of real stock returns which goes as /~°-^^,also shown. By using 
exponential weights (here of half-time corresponding to h ^ I0.5/T = 0.1), simulations 
yield a kurtosis that matches well to real stocks. 



Note that in our model, the behaviour of U (/) depends only trivially on 
the value of the parameter q (which essentially controls the strength of the 
feedback into the system) via the additive terms due to and the pre-factor 
bq which multiplies the i/i — yj terms. On a log-log plot, the slope of U{1) is 
constant for all q, the magnitude however vanishing as g — 1. If we normalize 
U{1) = 1 for / = for different values of q, then indeed all the U{1) curves 
collapse onto each other. Thus, it is not q but rather the temporal behaviour 
of the weights Wij that controls the rate of decay of U. In other words, the 
slow decay of volatility-volatility correlations is a signature of the memory of 
the system. 

Let us return our focus to the plot in Figure 2. We can compare the 
form of this volatility-volatility variogram with known results of real stock 
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'■■a Lag \n 

Figure 4: a) The average kurtosis decreases slowly across time scales following Eq(35). 
Here we keep T fixed and vary the relative half-time h = I0.5/T from h = 00 (equal weights 
with no decay) toward h = 0. We also show the line of slope —0.22 (a proxy for behaviom' 
of real stocks) and a line of slope —1 (a proxy for uncorrelated Gaussian variables) for 
comparison . (The fast decay as l/T increases is due to the cut-off, and is reduced if 
T ^ 00.) All curves are normalized to start at 1. b) This plot is a sub-set of the one 
above, together with two lines of slope —0.1 and —0.3. These lines quite well encompass 
the range of slopes of the curves for h = 0.02 — 0.1, shown also. On average, the slope is 
close to that observed for real stocks (~ —0.2). 

markets [1]. In particular, Bouchaud and Potters show that 

v={{dy^- dylif) <xv,- U{1)V2 (30) 

where vi and f 2 are case-specific constants. We aheady mentioned that they 
find that U{1) oc /"'^^ fits real stock returns quite well. This variogram was 
obtained from a simulation with /i = 0.1 which from the discussion in the 
above paragraph should reproduce a somewhat realistic decay of kurtosis, 
and therefore also a realistic variogram. Indeed, the solid line in Figure 2 is 
obtained from a fit of Eq(30) with coefficients vi = 1.8 and V2 = 1.73. The 
agreement is quite close. 

Figure 5 serves to illustrate yet again the feature of decreasing kurtosis. 
We have plotted out histograms of returns from a simulation of the process 
with q = 1.5 and an exponential half-time corresponding to h = 0.1. The 
distribution is tent-like for small time lags. As the time lag increases, the 
distributions become more Gaussian [1, 2, 26]. A systematic analysis of 
these distributions over time will be part of future work [29]. Finally, we 
also plot out the distribution of the instantaneous volatility calculated from 
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Figure 5: Shown here are histograms of returns yt+i — yt for different time lags I. From 
bottom to top corresponds the lags are I = 1,2, 4, 8, 16 and 32 multiples of the simulation 
time-step. For small / the distribution is highly non-Gaussian, becoming more Gaussian 
as the time scale increases. This is consistent with a decreasing kurtosis. 

the same simulation. As a proxy for the volatihty we have used squared 
returns averaged over 3 time increments. In Figure 6 we show the empirical 
distribution of the logarithm of the volatility from such a calculation. It is 
very close to a Gaussian distribution, as is quite apparent from the parabolic 
shape in this semi-log representation. This is encouraging, because it is 
well-known that real stock return volatility is well modeled by a log-normal 
distribution. 

It is worth mentioning that in principle one might be able to relax the 
assumptions of our initial model and still obtain processes with very similar 
features. Instead of arguing that at each time horizon i — j, the feedback 
process depends on Pg{yi \ yj), we could argue that it depends on Pg{yi \ yj) 
where yj can be any appropriate reference price, average of past prices, mean 
of the distribution of prices and so on. Exploring this in more detail is 
however beyond the scope of this paper. 
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Figure 6: The distribution of the effective volatihty appears similar to that of real stock 
returns, very close to a log-normal distribution, implying that the log-volatility follows a 
normal distribution. This is evident from the parabolic shape of the frequency plot shown 
here in semi-log representation. 

5 Conclusions 

In conclusion we have proposed a stochastic process with just one source 
of Brownian noise, yet the instantaneous volatihty is the result of a kind 
of statistical feedback across many time scales. This deterministic-volatility 
model captures many features observed in time series of real stock returns. 
These features or stylized facts include a log-normal distribution of volatility, 
together with volatility clustering and a kurtosis which decreases slowly over 
time. We find that the rate of this decay can easily be matched to that 
observed empirically for stock returns, and depends on the memory inherent 
in the system. 

Our model is parameterized by the index q, which corresponds to the 
degree of feedback into the system. For g = 1, there is no feedback and 
the process is identical to a standard Brownian motion. For g > 1 there is 
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a feedback, and there is a contribution to the instantaneous volatihty from 
multiple time horizons, such that the fluctuations depend on the probability 
of the latest observation at time i with respect to past observations on those 
different time horizons j. We model the probability densities as Tsallis (or 
Student) distributions, which they would be if each feedback process were 
independent. Mathematically, this non-Markovian process contains a sum 
over return differences weighted by a function of the time difference i — j- 
The exact form of this weighting controls the strength of the memory in the 
system. In this paper we explored the effects of both constant and expo- 
nentially decaying weights. For these two scenarios we could calculate the 
second moment as well as U{1), the decay of the fourth moment correlation. 
We found that on a log- log plot, the slope of this function U is constant 
for all q but varies according to the strength of the memory in the system. 
The function contains a multiplicative g-dependent pre-factor which goes to 
zero as q approaches 1, resulting in a vanishing kurtosis for that case. Other 
functional forms of the weighting kernel, such as a power-law or a weight 
which emphasizes natural time scales such as days, weeks or months, will be 
explored in future work [29]. 

These theoretical results in addition to the analysis of simulations showed 
that this model well-captures many stylized facts of real financial data. Fur- 
thermore it is simple and rather intuitive. The premise (that the volatility 
is a deterministic function of returns over different time horizons) is consis- 
tent with recent empirical analysis [27, 28]. In a forthcoming paper [29] , 
we shall analyze this and a related model in further detail, both theoreti- 
cally and with respect to empirical evidence of the underlying assumptions 
and possible predictions. A great challenge is to try and solve explicitly the 
distribution of the variables yt+i — yt- 

As a final remark we reiterate that in the case where there is only one 
relevant time, the current model reduces to the one which we previously 
proposed [20]. That model has been used for developing a theory of option 
pricing and another focus of future work will be to try and extend that 
theory to incorporate the multiple time horizons as introduced in the present 
framework. 
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A Appendix 

The quantity we would like to evaluate is Eq (29), namely 

U{l) = {^'){^li)-{Da^ + DC) (31) 

The cT^ terms are both equal to the diffusion coefficient At of Eq (24). 
We now replace the D, D and C terms with the sums they correspond to. 

i+l-N 

(^) = ( E <^-jV + K{^-j)-\yr-y3?) (32) 

j=i-N 

(D) = (£a,{z + l-jy + b,{z + l-j)-\y,+i-y,f) (33) 
j=i 

(C) ^ { f: a,{i + l-jr + b,{i + l-j)-\y,+i-yj)') (34) 

j=i+l-N 

We take the continuous limit of these sums, using the notation of Eq (21) To 
evaluate the integrals, we replace the {yi — yjY terms by their expectations 
according to Eq (24). We finally obtain 

U{1) = A| - (^^^ + (35) 



This result is valid for a constant equal weighted window. However it is 
straightforward to generalize for the case of exponential weights. Defining dyi 
as in Eq (20) we obtain (after taking the continuous limits and integrating): 

{D) = 6^AA(exp(-A(T-/)-exp(-AT)) 

+ ^(r(7 + i,A(r-0)-r(7 + i,A0) (36) 

(£>) = MA(l-exp(-AO) 

+ ^(r(7 + i)-r(7 + i,A/)) (37) 
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(C) = MA(exp(-A/)-cxp(-AT)) 

+ J(r(7 + i,A/)-r(7 + i,A(r-/)) (38) 

together with 

aX(r(- +i)-re,. + i.Ar)) 
(T^MA^ ^^^^ 

where r(a) is the standard Gamma function and r(a, x) represents the incom- 
plete Gamma function. Here r{a,x) corresponds to the incomplete Gamma 
function and r(a) the standard Gamma function. Inserting these equations 
into Eq(29) yields 

U{l)=Al-{D)A,-{D){C) (40) 

If we now let T go towards infinity in this calculation, we obtain results for 
the case of pure exponential decay. 
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